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1. CONSTRUCTING THE ANGULAR 
MOMENTUM 

We solve the bound state problem for two rel¬ 
ativistic fermions in the ladder approximation. 
The problem is studied in the framework of the 
explicitly covariant light-front dynamics In 

this approach, the state vector is defined on an 
hyperplane given by the invariant equation uj-x = 
0 with = 0. The standard light-front, reviewed 
in m, is recovered for oj = (1,0,0,—1). Because 
of explicit covariance, the wave functions for a 
given angular momentum are decomposed in a 
few covariant spin structures multiplied by scalar 
spin components. 

For two scalar particles with zero and non-zero 
angular momentum, this problem was solved in 
. The light-front deuteron wave function was 
calculated, in a perturbative way, in |^. Now 
we develop a general method, based on explic¬ 
itly covariant light-front dynamics, to solve the 
two-fermion bound state problem. We calcu¬ 
late the spin components non-perturbatively for 
one-boson exchange kernels: scalar, pseudoscalar, 
vector and pseudovector exchanges 

The interaction kernel U is a scalar, i.e., in 
c.m.-system, it depends on the scalar products 
of all available vectors, including the vector n = 
determining the light-front orientaion, that 
is: V = V(k,k',Si,S 2 ,n,M^). Here Si,S 2 are 
fermion spin operators. Thus, for the two-fermion 
system, the angular momentum operator com¬ 
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muting with the kernel reads: 

J = — X d/dk] — i[n x d/dn] -I- si -I- S 2 - (1) 

Besides, the three operators J^,Jz and = 

(n-J)^ commute not only with the kernel, but also 
with each other. Therefore, one should construct 
the eigenstates of J^,A^ (and J^): 

= J(J+1)^ (2) 

AV = (3) 

and solve the bound state problem for given 

(J, a). The equations for the functions with differ¬ 
ent (J, a) are not coupled with each other. Then 
for given J > 1, one should construct from the 
states with different a, a superposition satisfying 
the angular condition. 


1.1. The state with J = O’*' 

The basis for the state with J = O’*' contains 
two spin structures. The counting rule is based 
on the dimension of the spin matrix forming the 
two-fermion wave function: N = (2si -I- l)(2s2 + 
l)/2 = 2 X 2/2 = 2. The division by 2 takes 
into account the parity conservation. The vector 
ft enters in the wave function on the equal ground 
with the momenta of particles. In the reference 
system where ki + k 2 = 0, the wave function of 
two fermions in the state J = 0^ is represented 
as: 


= Vmwl^ip{k,fi)aywl^, 

\ ( fi r h 

ip{k, 1 


= A(h + '-^i^h 


V2 


sin0 


( 4 ) 
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where k = ki = —k^, k = k/k. The functions 
fi .2 depend on k and 9, where k-h = cos 6*. It is 
convenient to represent this wave function in the 
four-dimensional form, in order to use in calcula¬ 
tions the trace techniques of the Dirac matrices. 
Namely, 

<Pcr2Criikl,k2,p,UJT) = a/toUo- 2 (/^l) (5) 

with (p = fiSi + /2<5'2 and 


2^26 


^2 = 

Efc 

/ 2mtD 

^)75, (6) 

2^/2mk sin 9 

V ^-p 

where 

Lo = Uc 

= 7^7°, 

£k = 'Jk'^ + rn?. 


When ki + k 2 = 0, eqs. (i, (D turn into d)- 
The basis (^) is orthonormalized. 


1.2. The states with J = 1"*" 

The wave function of the J = I’*' state 
(deuteron, for example) is determined by six spin 
structures. They are split in two families with 
two spin structures for J = l+jO = 0 and with 
four spin structures for J = 1+, a = 1. 

One can easily check that the function ijPik^ n) 
satisfying the equation (^) with a = 0 is parallel 
to n, i.e., it satisfies the condition ijp = 

It has the following decomposition: 


n) 




92 


( 0 ) (T-{k cos 9 — n) 


sin0 


(7) 


The function ip^{k,n) corresponding to a = I 
is orthogonal to n, i.e., it satisfies the condition 
= 0. It is convenient to introduce the vec¬ 
tors orthogonal to ft: 


k± 


k — cos 9n 
sin0 


= a — {n-a)n. 


Then the function '0^ obtains the following form: 

tf^(.k,n) = 

+ (^^k±{k±-a±) - a±^ (8) 

+ 93^ \f^ i^-n) + 94^ )J^i[k X ft] 


Both wave functions d) and (||) can be also rep¬ 
resented in the four-dimensional form |^d], sim¬ 
ilarly to dl)- The functions /i ,2 for J = 0 and 
5 n 2 fo'' J = I'*’, a = 0 states are determined by 
a system of two equations. The functions 
for J = 1+, a = 1 state are determined by a sys¬ 
tem of four equations. The equations are given in 


2. SATISFYING THE ANGULAR GON- 
DITION 

The angular momentum operator in light-front 
dynamics is a dynamical one. Its equivalence to 
the kinematical operator (in particular, to the op¬ 
erator (|^ for two fermions) is ensured by the so 
called angular condition (see Q). This condi¬ 
tion cannot be satisfied for the state with fixed 
a 0, its solution is a superposition of the states 
with different a. Therefore, the wave function for 
J = I’*', satisfying the angular condition, is rep¬ 
resented as: 

tf{k,n) = coij^{k,n) + Ciijj^{k,n) (9) 

Since we don’t solve the full field theory, but re¬ 
strict ourselves by a particular (OBE) interaction, 
the angular condition can be also imposed in an 
approximate form. One can require that the wave 
function ip(k,ft) at A: = 0 does not depend on ft. 
The coefficients cq , ci in (||) are unambiguously 
determined by this requirement. 

The masses Ma for a = 0,1 are split, though 
they should be degenerate in the exact solu¬ 
tion. In standard LED, this split - the mass 
dependence on the angular momentum projec¬ 
tion - is the violation of the rotational invari¬ 
ance. Solution (j|) corresponds to average mass 
= CqMq -I- c\Mf. For two scalar particles 
with massless exchange, this formula gives results 
close to the Bethe-Salpeter solution and restores 
with high accuracy the Coulombien degeneration 

i- 

One can show that at fc = 0, the components 
of the wave function eq. (^, depend 

on 9 as: g^\k = 0,0) = bo cos 9, g^\k = 
0,9) = — 6osin0, where bo is a constant. In 
the wave function eq. (^, the component 
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g^i\k = 0, 0) = bi does not depend on 9, whereas 
g^li(k = 0, 0) = 0. The coefficients in the solu¬ 
tion (^) are expressed through bo^bi: 


Co 


bi 

y/2bl + bf 


Cl 


V^bo 

a /202 _|_ ^2 


( 10 ) 


The wave function is now can be rewritten as: 
if{k,n) = fi^a + f2^{3k{k-a) - a) 

+/3^(3n(n-<T) - <t) 

-|-/ 4 -( 3 fc(n-(T) -I- 3n{k-a) — 2{k-n)a) 

+f5i\l^[k xn]+ fo^{n{k-a) - k{n-a)). 


The six components fi are expressed through g^l 
and 5 ^/ 4 . 



Figure 1. The spin components fi{k,9 = 30°) 
determining the J = state wave function for 
scalar exchange. 

published in Q. 


3. NUMERICAL RESULTS 
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For illustration, we give here the results of nu¬ 
merical calculation with the kernel corresponding 
to the exchange of a scalar meson with mass g,. 
The form factor in the fermion-meson vertex is 
not introduced, but the calculation is done on the 
finite interval k < k^ax- Note that for J = O'*" 
and for the coupling constant a = g'^l{ATT) < 
Oc ~ 3.72 the solution is stable at kmax 00 
||. For J = 1+, m = 1, /I = 0.25 a = 1.18, 
B = 2m — M = 0.0501, kmax = 10, the en¬ 
ergy split between the states a = 0 and a = 1 is 
rather small: 2 %, in comparison to 20 % for two 
scalar particles |^. For two fermions interact¬ 
ing by pseudoscalar exchange, the energy split is 
much more significant: for a = 60, i?a=o = 0.103 
and Ba=i = 0.0494, that is more than a 100% 
splitting in energy j|]. 

The components //fc, 0 = 30°) are shown in 
figure ^ One can see that for scalar exchange 
the dominating component is fi (turning into the 
S-wave in the nonrelativistic limit). For the co¬ 
efficients the following values were found: 
Co = 0.5817, Cl = 0.8134, that is very close to 
1/V3 « 0.5773, ~ 0.8165. 

The details of approach and the numerical re¬ 
sults for all exchanges (S, PS, V and PV) will be 
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